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Abstract 

Topological properties of Fortuin-Kasteleyn clusters are studied on the 
torus. Namely, the probability that their topology yields a given subgroup 
of the first homology group of the torus is computed for Q = 1, 2, 3 and 4. 
The expressions generalize those obtained by Pinson for percolation (Q = 1). 
Numerical results are also presented for three tori of different moduli. They 
agree with the theoretical predictions for Q = 1, 2 and 3. For Q = 4 agreement 
is not ruled out but logarithmic corrections are probably present and they 
make it harder to decide. 

KEY WORDS: Potts models; bond percolation; torus. 

1 Introduction 

The Fortuin-Kasteleyn (FK) formulation of Potts models is often referred to as the 
Q-state bond-correlated percolation. The clusters formed by sites and bonds, called 
FK clusters, may wrap the torus in a non-trivial way. One can study the particular 
topology of a cluster by looking at its first homology group. This group will be 
isomorphic to a subgroup of the homology group of the surface on which the model 
is denned. In general, consider Q, the topological space given by all FK clusters of a 
configuration imbedded in a Riemann surface S. There is a natural homomorphism 
: Hi{Q) — > H\(S) between first homology group of each space. Therefore one 
may look for the probability that <p(Hi(Q)) be a given subgroup G of Hi(S). A 
definition of this observable was first given by Langlands et aim 0. In a preceding 
work, Pinson computed these probabilities in the percolation case on the torus using 
a method introduced by di Francesco et al §|. 

This paper generalizes Pinson's work on percolation to Potts models on the torus 
for Q — 1, 2, 3 and 4 to compute the probability that a FK configuration of a Q- 
state Potts model yields a particular subgroup of the first homology group of the 
torus. This observable is also studied numerically using Monte-Carlo simulations. 

The general concepts required for our calculation are defined in section 2. Then 
essential points on the transformation from the Potts model to a gaussian free field 
are discussed before we deduce the explicit expressions for probabilities in section 4. 
As expected, these are found to be modular invariant. Numerical results presented 
for three tori with different ratios are in good agreement with values computed 
except for Q — 4. For the latter case, it is difficult to estimate the continuum 
limit numerically because what looks like logarithmic finite-size corrections seem to 
appear. 
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2 Definitions 



A FK configuration (or FK subgraph) Q is built on a Potts configuration a by putting 
bonds between all same-state neighbors and then removing them with probability 
1 — p = 1 — exp(— K) where K is the critical ratio J/ksTc. The resulting graph Q 
is said a FK subgraph associated to a. Neither correspondences a — ► Q nor Q — > a 
are unique. A maximally connected component in is called a FK cluster. 

In this paper, we define the model on the torus, i.e. C/(Z + Zr), r G Z. The 
first homology group of the torus T is isomorphic to Z x Z. Therefore its subgroups 
are isomorphic to {0}, Z x Z and {m,n} = ((m, n)) where m, n G Z. 

Let a and /? be two independent cycles on the torus. On the parallelogram 0, 
r, 1, 1 + r, we choose a to be horizontal and (3 to be parallel to the edge (0, r). 
We distinguish different classes of FK clusters. They can be homotopic to a point. 
They may wind the torus a times along a and b times along (3, where a, b G Z and 
a A b =gcd(a, b) must be one. A cluster of this type will be called a cluster {a, b}. 
We define the signs of a and b as follows. Let us say we draw a boundary curve of 
such a cluster. Then a > if we go from left to right and b > if we go from top 
to bottom. Hence a cluster {a, b} is equivalent to a cluster {—a, — &}. The last class 
is formed by cluster with a cross topology, precisely clusters containing at least two 
connected subcomponents with different non-trivial windings {a, b}. 

We consider the natural homomorphism <p : Hi(Q) — > Hi(T) for a FK subgraph 
and the torus T. We are interested in the subgroup G of Z x Z which lie in <fi(Hi(Q)). 
Then, we associate to Q this particular subgroup G. In fact, the correspondence is: 

• G={0} if Q contains only clusters homotopic to a point; 

• G={a, b} if Q has only clusters homotopic to a point and at least one cluster 
of class {a, 6}; 

• G=Z x Z if contains a cross. 

Some examples are given in figure [I] for Q — 2. The two states are depicted re- 
spectively by black dots and white dots. Bonds were added between same-state 
neighbors with probability p. The subgraphs associated with each subgroup will 
be respectively called: trivial subgraphs, subgraphs {a, b} and subgraphs with a 
cross topology. Because of the nature of FK clusters the probability of G — {a, b}, 
a A b 7^ 1 is zero. Hence we will assume in the rest of the paper that a A b = 1. 

Let Zq(G) be the partition function restricted to FK configurations associated 
to the subgroup G in the continuum limit. The probability that a subgraph yields 
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ZxZ {0} {1,0} {1,-1} 




Figure 1: Different Q = 2 configurations with their associated subgroup. 
G is then 

„(G) = (1 ) 
where Z Q is Z Q ({0}) + Z Q (Z x Z) + £ aM=1 Z Q ({a, b}). 

3 Fortuin-Kasteleyn clusters and gaussian free fields 

To compute each weight, a tranformation of the Potts model in the FK formulation 
to the F model (or six- vertex model) and finally to the SOS model will be performed. 
The renormalization of the SOS model into a bosonic free fields is then used to get 
the probabilities in the continuum limit. 

The first step of the transformation involves a decomposition of subgraphs into 
oriented-polygon decompositions. There is a bijection between subgraphs and pos- 
sible non-oriented polygon decompositions on a lattice of¥xiV sites 0. The 
partition function of the Potts model in the FK formulation can then be rewritten 
as a sum of separated weights at the critical temperature ||. § : 

Z M,N = Z MXN {{Q}) + Z MXN {Z xZ)+ J2 Z^ N ({a, b}) 

aAb=l 

= Yl Q Np/2 +Q E Q Np/2 + E E ® Np/2 - ( 2 ) 

V{0} V{ZxZ} aAb=lV{a,b} 

where V{G} stands for all decompositions whose subgraph generates G. Here, iVp 
stands for the number of polygons in a specific decomposition. The Q factor in the 
restricted weight of Z x Z comes from a different Euler relation between the number 
of sites, polygons, connected components and faces for connected components with 
a cross topology. 
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One can build a bijective correspondence between trivial subgraphs and sub- 
graphs with a cross topology which preserves the number of polygons in the polygon 
decomposition. Let Q + be a subgraph of the cross topology type with decomposition 
V . A trivial subgraph Q is built on the dual lattice of Q + keeping the same V by 
putting its sites in the middle of each smallest square formed by four sites of Q+. 
Then, on the dual lattice (sites in the middle of plaquettes), a subgraph Q is con- 
structed by putting all bonds that do not intersect with bonds of Q + on the original 
lattice. The correspondence is bijective by construction and Go is necessarily trivial. 
One can also see that Q + and Qq have the same number of polygons. Consequently, 
we obtain an important duality relation for weights of {0} and Z x Z in (0): 

Z^ N (ZxZ)=QZ^ N ({0})- (3) 

In the F model, we express the factor Q Np / 2 in terms of a random orientation on 
each polygon. Let us write q = ±1 for non-trivial polygons (type n-h) and €j = ±1 
for polygons homotopic to a point (type h). We impose: 



Q N P /2 =Q N n _ h /2 Q N h /2 = J2 COS 




ex p 




{<=i=±l} 


i 


{6 3 =±1} 


j 



where eo and A are chosen such as Q 1 ^ 2 = 2cos(7re /2) and Q = 2cos(4A). We do 
not need the last details of the transformation to the F model for our purpose, but 
the reader could find them in 

To construct the SOS model from the F one, a field must be defined on the faces 
of the lattice || . These faces are formed by polygonal curves of type n-h or h. A 
polygonal line means an increase or a decrease by a factor 7r/2 of the field depending 
on the polygon orientation. Discontinuities or frustrations are therefore induced in 
the field by polygons n-h of the decomposition. For subgraphs with a cross topology 
and trivial one, only trivial polygons are present; thus the field is always periodic. 
For subgraphs {a, b}, if we write respectively 8(j>i and 5<p T for discontinuities observed 
along the cycles a and /3, we have : 

50i = i^ )S ^ J e i = 11111 <50r = "^ a ^^ e « = 7rm ' (5) 

i i 

where m, m! G Z and the sum is only over n-h polygons. For subgraphs wih a cross 
topology and trivial one, m and m' are zero. Subgraphs {a, b} have possibly polygons 
h and necessarily two polygons n-h of class ±[o, b] for each non-trivial connected 
component. Therefore, their boundary conditions will be: m = bk, m' = ak for a 
certain fceZ. Moreover 
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Remark 1 because aAb = 1, only a unique subgroup {a, b} can produce non-periodic 
field with given discontinuities m, m' , (m, m') 7^ 0. 

We can express the correction term due to non-trivial polygons in (|J) as cos [|e J2i e «] = 
cos[7re (m A m')] with the help of (||) and the fact that a Ab = 1. 

The extra factor Q in front of the second term in (|2|) does not appear in the usual 
partition function of the F model. We have not taken it into account yet in the SOS 
model partition function. Consequently, let us first consider Zq IxN , the partition 
function of the Q-state Potts model with a weight Q times lower than the correct 
one for subgraphs with a cross topology. In regard to the transformation performed, 
it corresponds to the SOS model partition function with the non-trivial correction. 
Considering that the SOS model is assumed to renormalize onto a gaussian free field, 
Zq IxN becomes in the continuum limit ||: 



Zq xN (eo) - Z(g, eo) = ^ Z m>m ,(g/4) cos [7re (m A m')) 



(6) 



m,m'i 



The hat is the reminder of the missing factor Q. The value of g is given by Q = 
2 + 2 cos (irg/2), g G [2,4] (see 0). The term Z m m i{g) is the bosonic partition 
function on a torus r = r^+ir/ for a field with discontinuity 81(f) = 2irm, 5 T <p = 2nm' 
and coupling constant g: 



V9 



exp 



-irg- 



m 2 rf + (m' — ruTR)' 



(7) 



where q = exp(2i7rr) and 77(g) is the Dedekind function. 



4 Probabilities of homology subgroups 

Each specific weight in (0) can be derived using expressions of the last section. The 
partition function is then easily computed to get the probability ([I]) for a particular 
subgroup G. Every expression found must be modular invariant. 

We first find the weights of {0} and ZxZ. We notice that for eo = 1 in (^), the 

quantity Q Nn - h l 2 is zero , therefore Zq({o, 6}) 

we get: 



= 0. Hence, if we use this in 

e =l 



%l)=Z o ({0}) + ^(ZxZ). (8) 
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The desired results then follow from (^) 

^q({0}) = - ^ Z mim /(p/4)cos[7r(m Am')] 



m,m'£X 



Z Q (ZxZ) = Q ^ Z m,m'(g/4) cos [n(m Am')]. 



(9) 
(10) 



m,m'£l 



We introduce the notation Z^' (g/A) for the weight of all subgraphs {a, b} which 
induce periodic field and write Z mtm i(g/A,e ) = Z mim /(g/A) cos [7re (m A m')\ for 
convenience. The weight of the subgroup {a, b} is then given by the sum over all 
non-zero boundary conditions particular to this subgroup by remark |I| plus the 
contribution of subgraphs {a, b} with periodic field: 



11) 



m=bk 
m'=ak 
k£2 



For eo = 1, it follows directly from (|TT|) that: 



(12) 



m=bk 
m'=ak 

kez\{o} 



We insert this new expression for Z^'^g/A) in (|TT). The correct weight is then given 
by: 



z Q({a, b}) = ^2 [ z m,m'(g/^ eo) - Z m>ml (g/A, 1) 



(13) 



fcez\{o} 

The partition function is computed by summing over all subgroups: 

z q = y^-z( g ,i)+ E [ix(9/ 4 . e o)-iy(s/4,i: 

aAfe=l m=bk 
m'=ak 
k€Z\{0} 

= Z(g,e ) + ^Z^-Z[g,l] 



(14) 



in agreement with reference ||. We added zero to the r.h.s of the first equality in 
([14]) by considering the term m = m' = in the sum. 
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Each probability can then be calculated with ([I]). For a specific Q, sums over 
Z(g/A, e ) involved in the expressions found are computed using restricted sums over 
equivalence class of m A w! so the cosinus is a constant in each of them. 

The subgroups {0} and Z x Z are both stable under the action of any modular 
transformation of the torus whose generators are T : r — ► r + 1 and S : r — > 
— 1/r. Hence, their probability must be invariant. For subgraphs {a, b}, modular 
transformations of non-trivial curves impose: 



tq({M)}) 
tq({M)}) 



n Q ({(a + b, &)}) 
n Q ({(-b,a)}) 



T+l 



-1/7 



(15) 
(16) 



These restrictions are easily verified directly on ([13]) with the help of the modular 
transformation of the bosonic partition function H: 



J ra,m' 



J m,m+m' 



Zm'—m 



T + l 



-1/r 



(17) 
(18) 



In the particular cases Q — 1 and Q = 2, the expressions for 7TQ({a, b}) are 
written in an elegant way using theta functions H . The calculations yield for r a ^ = 
a — 6t 



7n({a,&}) 
7r 2 ({a,6}) 



1 



2|r a , b ||r / (g)| 2 
1 

\Ta,b\\v((l)\ 



9, 



IT! 



6 l r a,fc| 2 



0, 



3iri 
2\r n J 2 



29o 



3iri 
2|r a , b | 2 



2T/ 

2 1 3|r^ 



73 1 sCT 



2 (r) + e 3 (r) + 9Jt) 



(19) 
(20) 



5 Numerical simulations 

Numerical simulations were performed on square lattices with V x H sites on three 
different tori, r = i, r = 2i and r = i + 1/2. We chose to consider the quantity 
7Tq({0}) + 7Tq(Z x Z) instead of two different values because they are simply related 
by (|3p for every lattice. This allowed us technically to experiment on larger samples. 
Results are shown in Table [I]. 

On a V x H lattice, the identification rules on the boundaries for a site 
were (i, H ) = (i,0), (V,j) = (0,j) for r = % and r = 2i and (z, if) = (i,0), 
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(y,j) = (0, j + V/2 mod V) for r = i + 1/2. We used Swendsen-Wang algorithm 
for Q = 2, 3 and 4 with at least 8000 initial thermalisations and five upgrades 
between measurements. 

We proceed differently to approximate the continuum limit for Q — 1, 2 than 
for Q — 3, 4. In the first case, we took lattices of dimensions 256 x 256 sites for 
t = i and r = i + 1/2, 256 x 128 for r = 2i. These are indeed good approximations 
of the continuum limit as can be seen by comparing measurements on lattices with 
half this linear size. In Table 1, the numbers after the vertical bar represent the 
statistical error on the last two digits for a 95%-confidence interval. For example, 
0.0210|10 means the interval [0.0200,0.220]. Each probability was calculated on 
a sample larger than 2 x 10 7 configurations. The figures in bold faces are the 
predictions obtained from (Q), (^TJ) and (|TB|). The agreement between numerical 
data and analytic calculations are convincing. 

For Q = 3 and Q = 4, a more complex method is necessary because the scale 
effect is still visible in large lattices. Consequently, we used a power law hypothesis 
& 

\n Q (G,H,V)-7t Q (G)\~aVP. (21) 

where a > and (3 < are two parameters to be determined. We obtained 
a and (3 by a linear regression of the function f[i] = log \71q(cl, b, 2 l H , 2 l V ) — 
7r Q ({(a, 6)}, 2 i+1 H , 2 i+1 V )|. We took V = 16 for Q = 3 and V = 32 for Q = 4 with 
i = 0,1, 2, 3, 4. The errors were estimated using a numerical method. For each i, we 
had a sample of at least 2 x 10 7 configurations. The results are very good for Q = 3 
although errors are larger using this method. One might ask whether the exponent 
(3 of ( |2~1"D is universal. However, the values of the exponent (3, say in the case Q = 3, 
differed for the different subgroups and the different ratios studied though they re- 
mained of order 1. For Q = 4, we naively tried the same approximation. These 
results in Table 1 and the behaviour plotted in figure |2| suggest that the finite-size 
correction term is not of the form (F2T|). Several authors have argued convincingly 
that finite-size effects for the Q = 4 Potts model show logarithmic corrections [|J. 
It is hard to see them in our data since the fit would involve more parameters and 
thus would need more points to be accurate. Thus one would have to experiment 
on lattices larger than 512 x 512 and to deal with the problem of getting samples 
with a relatively small statistical error. 
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Figure 2: Deviation of numerical results from theoretical prediction for 7r(Z x Z) + 
7r(0) on a N x N lattice (r = i) as a function of log 2 (A r ) 
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Table 1: Probabilities calculated with both analytical and numerical method. 



Q 


T 


7r Q ({0}) + 7r o (ZxZ) 


t«({1,0}) 


ttq({0,1}) 


t«({M» 


ttq({1,-1}) 


1 


i 


.61908 


.16942 


.16942 


.020980 


.020980 






.61903|22 


.16948(17 


.16939(17 


.020995(64 


.020957| 64 




2i 


.45376 


.50304 


.024950 


.008755 


.008755 






.45380|22 


.50293(22 


.024983(70 


.008749(42 


.008779(42 




i + 1/2 


.62840 


.16816 


.10006 


.10006 


.001519 






.62849|22 


.16812(17 


.10008(13 


.09998(13 


.001504(17 


2 


i 


.67927 


.14644 


.14644 


.013903 


.013903 






.67907|21 


.14654(16 


.14645(16 


.013949(53 


.013940(53 




2% 


.47863 


.49625 


.015376 


.004734 


.004734 






.47888(32 


.49595(32 


.015384(78 


.004751(43 


.004760(43 




i + 1/2 


.69083 


.14569 


.08095 


.08095 


.000727 






.69091(22 


.14569(17 


.08100(13 


.08081(13 


.000724(13 


3 


i 


.74533 


.11867 


.11867 


.008660 


.008660 






.7457|21 


.1173(19 


.1187(14 


.00863(42 


.00871(49 




2i 


.52571 


.45995 


.0092445 


.0024968 


.0024968 






.5250|21 


.4605(21 


.00921(45 


.00240(29 


.00247(23 




i + 1/2 


.75821 


.11820 


.061448 


.061448 


.00033 






.7584(19 


.1180(15 


.0606(13 


.0608(13 


.00000(81 


4 


i 


.85034 


.071728 


.071728 


.0030997 


.0030997 






.8265(23 


.08343(16 


.08151(16 


.00380(43 


.00440(32 




2i 


.63368 


.36172 


.0032494 


.00067549 


.00067549 






.6013(25 


.3929(25 


.00369(28 


.00087(16 


.00085(14 




i + 1/2 


.86316 


.07151 


.03260 


.03260 


.000061 






.8426(23 


.0794(19 


.0389(11 


.0383(13 


.00011(4 
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